## irr_finding


from PyM import *

K = Zn(2)

[ZX,X] = polynomial_ring(K,'X')

show( is_irreducible(X**3+X**2+1))

show( is_irreducible(X**3+X+1))

f = sum(geometric_series(X,7))

show( is_irreducible(f))

show( factor(f))

def irred(r): return get_irreducible_polynomial(K,r)

show(irred(4))